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Abstract.  A  case  study  is  presented  of  6.5  h  of  simulta¬ 
neous  colocated  stratospheric  Differential  Absorption  Lidar 
(DIAL)  measurements  of  ozone  concentration  and  Rayleigh- 
Mie  Doppler  (RD)  lidar  measurements  of  horizontal  wind 
velocity  front  the  Observatoire  de  Haute  Provence  (OHP), 
France  (44°N,  6°E).  The  RD  lidar  observations  reveal  a  dis¬ 
tinct  gravity  wave  motion  at  ~  12-1 8  km,  while  the  DIAL 
data  show  laminated  ozone  structures  in  the  same  height 
range.  We  combine  these  data  to  show  that  advection  due  to 
the  gravity  wave  cannot  produce  the  large  ozone  lamina  at 
~14km. 

Introduction 

Stratospheric  ozone  abundances  exhibit  appreciable  meso¬ 
scale  variability.  Large  perturbations  with  vertical  scales 
~  1-5  km  are  observed  routinely  in  ozonesonde  data  [Ehhalt 
etal,  1983 ;  Hofmann  et  al,  1989;  Reid  and  Vaughan,  1991; 
Teitelbaum  et  al.,  1994;  Mitchell  et  al.,  1996],  while  mea¬ 
surements  from  stratospheric  aircraft  show  variability  over 
a  broad  range  of  horizontal  scales  [Nastrom  et  al..  1986; 
Danielsen  etal.,  1991;  Bacmeister  et  al.,  1996].  Since  ozone 
in  the  lower  stratosphere  (below  ~25  km)  is  photochemically 
inactive  on  timescales  up  to  5  days  (in  the  absence  of  hetero¬ 
geneous  processing),  then  such  variability  must  in  general 
have  a  dynamical  origin. 

Gravity-wave-induced  vertical  advection  was  first  sug¬ 
gested  as  a  possible  cause  of  ozone  perturbations  by  Chiu  and 
Ching  [  1978] .  Aircraft  observations  during  the  Stratosphere- 
Troposphere  Exchange  Project  (STEP)  revealed  inertia-grav¬ 
ity  waves  (IGWs)  and  seemingly  related  perturbations  of 
ozone  and  other  constituents  [Chan  etal.,  1991;  WilsonetaL, 
1991;  Alexander  and  Pfister,  1995].  Danielsen  et  al.  [1991] 
argued  that  the  vertical  and  horizontal  displacements  of  these 
IGWs  produced  the  observed  ozone  perturbations,  and  that 
irreversible  constituent  transport  resulted  as  these  waves  dis¬ 
sipated.  Teitelbaum  et  al.  [1994],  Alexander  and  Pfister 
[1995]  and  Teitelbaum  et  al.  [1996]  argued  that  gravity-wave- 
induced  vertical  advection  accounted  for  the  ozone  fluctua¬ 
tions  in  their  data. 
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Conversely,  Ehhalt  et  al.  [  1983]  found  that  long-term  vari¬ 
ances  of  ozone  inferred  from  ozonesonde  releases  at  Hohen- 
peissenberg  (47°  N,  1 1°E)  were  significantly  larger  than  those 
expected  on  the  basis  of  wave-induced  potential  temperature 
fluctuations  alone,  indicating  that  gravity  wave  vertical  ad¬ 
vection  cannot  explain  all  of  the  structure.  Later  analysis  of 
variability  in  a  large  base  of  ozonesonde  data  from  middle 
to  high  northern  latitudes  by  Reid  et  al.  [1993]  identified 
the  largest  perturbations  (so-called  “laminae”)  with  differen¬ 
tial  transport  across  the  polar  vortex  edge.  High-resolution 
parcel  advection  models  of  this  transport  have  succeeded  in 
simulating  laminated  vertical  ozone  profiles  which  are  sim¬ 
ilar  to  those  observed  [e.g.,  Orsolini  et  al.,  1995],  Indeed, 
Newman  and  Schoeberl  [  1 995]  recently  used  such  a  model  to 
argue  that  quasi-horizontal  flow  in  and  out  of  the  polar  vor¬ 
tex  could  explain  those  laminar  ozone  structures  observed  in 
STEP  data  which  Danielsen  et  al.  [1991]  had  previously  at¬ 
tributed  to  IGW  motions.  Recent  analysis  of  a  large  amount 
of  airborne  stratospheric  ozone  data  by  Bacmeister  et  al. 
[1996]  concluded  that  neither  IGWs  nor  large-scale  advec¬ 
tion  alone  could  explain  the  mean  mesoscale  variability  ev¬ 
ident  in  these  data,  and  that  both  effects  may  combine  and 
interact  to  produce  the  observed  variability. 

Recent  studies  by  Reid  et  al.  [1994]  and  Langford  et  al. 
[  1 996]  combined  velocity  and  ozone  data  to  investigate  IGW 
activity  in  ozone.  However  the  poor  temporal  resolution 
of  sonde  data  was  a  limiting  factor  in  the  study  of  Reid 
et  al.  [1994]  and  the  spatial  separation  of  75  km  of  the  radar 
and  DIAL  instrument  in  the  experiments  of  Langford  et  al. 
[1996]  meant  that  waves  observed  in  one  data  set  may  not 
always  relate  to  those  in  the  other.  Here,  we  report  on  a 
pilot  experiment  which  was  devised  to  avoid  many  of  these 
difficulties.  Two  colocated  lidar  instruments  were  used  to 
acquire  simultaneous  stratospheric  velocity  and  ozone  data 
with  high  vertical  and  temporal  resolution.  We  use  these 
velocity  data  to  identify  and  characterize  an  IGW  in  the 
stratosphere  and  to  predict  its  effect  on  ozone  profiles,  then 
compare  these  predictions  with  the  contemporaneous  ozone 
data  to  infer  the  processes  responsible  for  the  structure  in  the 
ozone  data. 


Observations 

The  DIAL  ozone  and  RD  lidars,  located  approximately 
10  m  apart  within  the  same  building  at  OHP,  were  run  si¬ 
multaneously  during  the  observation  period  (2230  LT  on  23 
October  to  0500  LT  on  24  October,  1995).  Table  1  sum- 


1627 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1997  to  00-00-1997 

4.  TITLE  AND  SUBTITLE 

Dual  lidar  observations  of  mesoscale  fluctuations  of  ozone  and  horizontal 
winds 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,E.O.  Hulburt  Center  for  Space 

Research, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

see  report 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

4 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


1628 


GIBSON- WILDE  ET  AL.:  MESOSCALE  OZONE  FLUCTUATIONS 


Table  1.  DIAL  and  RD  lidar  capabilities  at  OHP. 


Ozone  DIAL 

RD  lidar 

Beam 

vertical 

vertical  +  off-zenith 

direction 

40°  north  and  east 

Wavelengths 

absorbed  308  nm 
reference  355  nm 

532  nm 

PRF 

50  Hz 

30  Hz 

Height  resol. 

600  m 

115m 

Time  resol. 

30  min 

30  min 

Height  range 

12 -30  km 

12-45  km 

Major  error 

water  vapour 

aerosols 

sources 

aerosols 

inarizes  the  operational  specifics  of  both  lidar  systems  (see 
also  Gamier  and  Chanin  [1992]  and  Lacoste  et  al.  [1992]). 
Measurements  of  lower  stratospheric  (12-30  km)  ozone  con¬ 
centration  were  provided  by  the  DIAL  system  and  zonal  and 
meridional  velocities  were  provided  by  the  RD  lidar. 

Figure  1  presents  the  sequence  of  DIAL  ozone  concentra¬ 
tion  profiles  in  two  consecutive  plots.  The  estimated  error  in 
measurement  is  ~  10%  in  the  height  range  12-20  km.  The 
RD  lidar  zonal  and  meridional  velocity  profiles  in  Figure  2 
are  displayed  at  a  vertical  grid  spacing  of  150  m,  obtained 
from  regridding  using  spline  interpolation  from  the  initial 
115  m  resolution,  for  compatibility  with  the  DIAL  profiles. 
The  RD  lidar  data  are  acquired  typically  with  a  time  resolu¬ 
tion  of  5  min  so  that  the  30  min  profiles  presented  here  have 
an  estimated  error  of  less  than  5%  between  1 2-20  km. 

Figure  3  displays  the  average  ozone  concentration  profile 
calculated  over  the  6.5  h  observation  period.  The  error  asso¬ 
ciated  with  this  average  profile  is  ~2%  in  the  height  range 
12-20  km.  The  small-scale  features  evident  in  this  average 
profile  must  result  from  dynamical  activity  on  timescales 
greater  than  6.5  hours. 

Ozone  Laminae 


MERIDIONAL  WIND  PERTURBATIONS  (m  s  ’) 

Figure  2.  RD  lidar  profiles  from  2230  LT  on  23  October  to 
0500  LT  on  24  October.  1995.  The  zero  point  of  each  profile 
is  labeled  with  the  starting  time  of  its  30  min  observational 
integration.  Successive  profiles  are  offset  by  10  m  s-1. 


2.5  km  in  vertical  extent.  Reid  and  Vaughan  [1991]  found  that 
the  majority  of  laminae  in  mid-latitude,  northern  hemispheric 
ozonesonde  data  were  within  the  height  range  12-18  km. 

The  6.5  h  mean  lidar  ozone  concentration  profile  (Figure 
3)  was  compared  with  an  ozone  concentration  profile  from  an 
ozonesonde  launched  from  OHP  at  1200  LT  on  24  October 
1995,  7  h  after  the  final  lidar  profile  (not  shown  here).  The 
ozonesonde  profile  confirmed  the  lidar  observation  of  the 
lamina  at  14  km  so  that  a  minimum  lifetime  of  14  h  was 
estimated  for  the  lamina  at  14  km  based  on  its  occurrence  in 
both  the  lidar  and  ozonesonde  profiles.  The  lamina  at  18  km 
was  reduced  by  the  time  of  the  ozonesonde  flight.  Thus,  from 
the  lidar  observations  alone,  the  estimated  minimum  lifetime 
for  the  18  km  lamina  was  6.5  h. 


The  two  local  maxima  in  ozone  concentration  in  Figure 
3  at  ~14  km  and  ~18  km  are  examples  of  ozone  laminae. 
They  were  identified  as  such  using  the  definition  of  Reid  and 
Vaughan  [1991]  which  requires  large  positive  or  negative 
ozone  concentrations  which  are  limited  to  a  layer  less  than 


0200  0300  0400  0500 

OZONE  CONCENTRATION  (1012  cm  3) 


Figure  1.  DIAL  ozone  concentration  profiles  from  2230 
LT  on  23  October  to  0500  LT  on  24  October,  1995.  Each 
profile  is  labeled  by  the  start  time  of  its  30  min  integration. 
Successive  profiles  are  offset  by  5x  1012cm-3. 


An  Inertia  Gravity  Wave 

Successive  vertical  profiles  of  zonal  and  meridional  ve¬ 
locities  from  the  RD  lidar  (Figure  2)  indicate  a  wave-like 
feature  below  20  km  with  downward  phase  (upward  energy) 
propagation  which  persists  throughout  the  6.5  h  observation 
period.  A  hodograph  analysis  (Figure  4)  shows  mostly  clock¬ 
wise  rotation  of  the  velocity  fluctuations  with  height  and  a 
well-defined  velocity  ellipse,  which  suggests  that  this  feature 
is  a  long-period  upward-propagating  IGW.  The  wave  param¬ 
eters  were  derived  from  these  data  using  a  Stokes  parameter 
analysis  [e.g.,  Eckermann,  1996]  and  are  presented  in  Table 
2.  The  values  of  observed  and  intrinsic  wave  frequency  were 
used  to  infer  that  the  horizontal  propagation  direction  of  the 
IGW  was  towards  the  south-west. 

Vertical  IGW  Displacements 

The  IGW  parameters  in  Table  2  were  used  to  estimate  the 
influence  of  vertical  advection  due  to  this  IGW  on  the  ozone 
concentration  profile.  For  adiabatic  motion  and  a  constant 
background  ozone  gradient,  the  relative  ozone  concentration 
perturbation,  O3/O3,  is  given  by  [Chiu  and  Ching,  1978] 


<%(*.*)  „  (d03(z) 

Oi(z)  ~  ^  Oi(z)  )\  dz  7 H  )' 


(1) 
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Figure  3.  Mean  6.5  h  profile  of  ozone  concentration  (bold 
line)  and  fitted  mean  profile  (dashed-dotted  line).  A  5%  error 
estimate  is  indicated  by  the  dotted  line. 


where 

C'(*,i)«2^M«i/2(i>)  (2) 

is  the  vertical  displacement  of  a  parcel  by  the  IGW,  u  is  in¬ 
trinsic  frequency,  <5_  (a>)  —  1  —  (//cD)2,  /  is  the  inertial  fre¬ 
quency,  U'(z,t)  is  the  horizontal  velocity  perturbation  along 
the  wave  vector,  030)  is  the  background  ozone  concentra¬ 
tion  (assumed  time-invariant),  H  is  the  pressure  scale  height, 
7  is  the  ratio  of  specific  heats  (~1.4),  and  N  is  the  Brunt- 
Vaisala  frequency  (~0.02  rad  s-  1  based  on  ozonesonde  tem¬ 
perature  data).  The  030)  profile  was  estimated  by  a  low- 
pass  filtered  cubic  polynomial  fit  to  the  mean  profile,  which 
removes  the  laminae  (see  Figure  3). 

Figure  5a  shows  the  measured  03(z,  t)/03(z)  deviations 
from  the  cubic  polynomial-fitted  profile,  as  well  as  the 
expected  perturbation  produced  by  the  IGW  in  Figure  4, 
calculated  by  first  converting  U'(z,t)  to  a  ( '(z,t)  oscilla¬ 
tion  using_(2)  (see  Figure  5b),  then  using  this  to  calculate 
030>O/030)  using  (1).  The  relative  magnitudes  of  the 
measured  and  calculated  ozone  perturbations  in  Figure  5a 
show  that  the  vertical  advection  produced  by  the  IGW  in 
Figure  4  can  drive  only  a  5%  relative  ozone  fluctuation,  ex¬ 
cept  in  the  height  region  below  ~  14.5  km  where  the  inverse 
dependence  on  030)  in  (1)  leads  to  perturbations  of  up  to 
10%.  The  estimation  of  030)  using  the  smoothed  profile 
with  laminae  removed  is  the  key  factor  in  this  calculation, 


Figure  4.  Hodograph  between  12-20  km  of  mean  6.5  h  zonal 
and  meridional  wind  fluctuations  smoothed  over  600  m.  The 
bold  line  is  the  major  axis  of  the  gravity-wave  ellipse  from 
the  Stokes  parameters  analysis. 


Table  2.  IGW  parameters  derived  from  successive  RD 
lidar  profiles  and  Stokes  parameters  analysis. 


Wave  Parameter 

Estimated  Value 

Vertical  wavelength 

2.3  km 

Observed  vertical  phase  speed 

-0.06  m  s-1 

Observed  frequency 

1.6/  ~  0.00016 rads-1 

Intrinsic  frequency 

2.2/ ~  0.00022  rads-1 

Horizontal  wavelength 

360  km 

and  gives  the  expected  response  in  the  absence  of  preexist¬ 
ing  ozone  lamination.  Figure  5  clearly  illustrates  that  the 
lamina  at  14  km  cannot  result  from  vertical  advection  by  the 
coexisting  IGW  in  Figure  4. 

Horizontal  IGW  Displacements 

Horizontal  IGW  displacements  were  neglected  in  (1)  be¬ 
cause  the  lidar  data  do  not  provide  information  on  mean 
horizontal  gradients  in  03.  However,  when  zonal  and  merid¬ 
ional  03  gradients  are  nonzero,  additional  contributions  to  (1) 
arise  due  to  the  zonal  (%')  and  meridional  (if1')  displacements 
of  the  wave  [e.g.,  Danielsen  et  al.,  1991],  Since  the  IGW 
propagates  south  westward,  with  u  «2.2/  and  W  peaking  at 
~4m  s-1,  then  x'  and  if/  both  oscillate  air  parcels  ±12.7  km 
zonally  and  meridionally  from  their  equilibrium  positions. 

For  horizontal  IGW  displacements  to  produce  the  lamina 
at  14  km,  large  increases  in  03  must  occur  no  more  than 
~20  km  away  from  the  lidar  site  so  that  the  IGW  can  advect 
these  larger  ozone  abundances  into  (and  then  out  of)  the  lidar 
field-of-view.  This  also  means  that  such  an  IGW-produced 
lamina  could  exist  over  a  given  location  for  no  longer  than 
one  ground-based  period  of  the  wave,  which  from  Table  2 
is  ~10.S  h.  Since  the  lamina  at  14  km  was  observed  to  be 
stable  over  at  least  14  h,  it  is  clear  that  horizontal  advection 
due  to  the  IGW  cannot  explain  this  feature. 

Since  advection  by  the  IGW  cannot  produce  the  observed 
lamina  at  14  km,  this  strongly  suggests  that  it  resulted  from 
synoptic-scale  advection  patterns.  Initial  indications  from 
isentropic  parcel  advection  calculations  suggest  that  a  fila¬ 
mentary  structure  was  indeed  present  at  approximately  14  km 
over  OHP  during  this  experiment  [R.J.  Atkinson ,  private  com¬ 
munication].  The  lamina  at  18  km  appeared  to  vary  in  shape 
and  intensity  from  profile  to  profile  (see  Figure  1).  Con¬ 
sequently,  IGW-induced  horizontal  advection  could  play  a 
more  significant  role  in  driving  this  feature. 

Conclusions 

The  dual  lidar  measurement  technique  described  here  per¬ 
mits  a  more  thorough  study  of  mesoscale  dynamics  and  tracer 
variability  than  has  been  possible  previously.  The  RD  lidar 
data  revealed  an  upward-propagating  IGW  below  20  km, 
with  a  vertical  wavelength  of  2.3  km,  an  observed  vertical 
phase  speed  ~0.06  m  s'"1,  and  a  peak  horizontal  velocity 
amplitude  of~4ms_1.  We  were  able  to  show  conclusively 
that  the  strength  and  persistence  of  an  ozone  lamina  at  14  km 
in  DIAL  data  was  not  explained  by  wave  advection.  Con¬ 
versely,  the  observed  variability  of  a  lamina  at  18  km  was 
somewhat  more  consistent  with  IGW  influences. 

The  time-height  coverage  and  resolution  of  both  data  sets 
were  vital  in  allowing  us  to  form  these  conclusions.  Further 
observations  using  this  technique  should  be  especially  valu¬ 
able  in  investigating  gravity  waves,  filamentary  ozone  struc¬ 
tures,  and  their  interactions,  which  Bacmeister  et  al.  [1996] 


1630 


GIBSON- WILDE  ET  AL.:  MESOSCALE  OZONE  FLUCTUATIONS 


-604020  0  204060  -300200100  0  100200 

RELATIVE  OZONE  (%)  VERT.  DISPLACE,  (m) 


Figure  5.  Measured  relative  ozone  concentration  pertur¬ 
bation  determined  by  subtracting  a  cubic  polynomial  mean 
profile  (dotted-dashed  curve  in  panel  (a))  and  calculated  re¬ 
sponse  of  the  IGW  using  (1)  and  (2)  (solid  curve  in  panels 
(a)  and  (b),  respectively). 

have  argued  may  be  very  important  in  explaining  observed 
power  spectra  of  ozone  variability  in  the  lower  stratosphere. 
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